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ABSTRACT 

New astrometric analysis of four WFPC2 images of the isolated neutron star 
RX J185635-3754 show that its distance is 117±12 pc, nearly double the originally 
published distance. At the revised distance, the star's age is 5x10^ years, its 
space velocity is about 185 km s~^, and its radiation radius inferred from thermal 
emission is Roo ~ 15 km, in the range of many equations of state both with and 
without exotic matter. These measurements remove observational support for an 
extremely soft equation of state. The star's birthplace is still likely to be in the 
Upper Sco association, but a connection with C, Oph is now unlikely. 

Subject headings: stars: neutron; stars: individual (RX J185635-3754); stars: 
kinematics; open clusters and associations: Sco-Cen 

1. Introduction 

The compact object RX J185635-3754 (Walter, Wolk, & Neuhauser 1996; Walter & 
Matthews 1997) is one of the closest isolated neutron stars to the Sun (Walter 2001). Be- 
cause of its relative brightness, its isolated nature, and its apparently thermal spectrum 
(Pons et al. 2002) from X-ray to optical wavelengths, this object affords the opportunity 
to study the surface emission properties of neutron stars, and to measure its radius. These 
are important constraints on the dense matter equation of state and the interior compo- 
sition of the neutron star. The angular diameter has been estimated by modeling of the 
spectral energy distribution (e.g.. Pons et al. 2002), while the distance is inferred from the 
trigonometric parallax. 

Walter (2001) found a parallax of 16.5±2.3 milli-arcseconds (mas), based on three images 
obtained with the Hubble Space Telescope Wide Field and Planetary Camera 2 (WFPC2). 
The implications of the inferred 60 pc distance include a likely origin in the Upper Scorpius 
OB association, possibly as a companion to the runaway O star C, Oph, an age of 0.9x10^ 
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years, and a space velocity of about 100 km s ^. The uniform temperature heavy-element 
atmospheric models of (Pons ct al. 2002) yield a radiation radius = R/ \/l — 2GM/c^ ^ 
6 - 8 km and a redshift z = {1 - 2GM / Rc^)'^/'^ - 1 «i 0.3 - 0.5, where M and R are the 
neutron star's mass and radius. The indicated radii and masses are in the ranges 4.5-8 km 
and 0.6-1.2 M0, respectively. The implied radius R is smaller than allowed by any reasonable 
equation of state, including that of self-bound quark matter, unless the mass is less than 0.8 
M0. In that case the radius could be matched by self-bound quark matter configurations. 
Relaxing the assumption of uniform surface temperature. Pons et al. (2002) showed that 
somewhat larger radii, up to Roo=13 km, are accommodated. The inferred small radius, 
together with the lack of photospheric features in the X-ray spectrum (Burwitz et al. 2001) 
and the lack of strong pulsations (Ransom, Gaenslcr, & Slane 2002), has led to speculation 
that the object might be a self-bound quark star (e.g., Xu 2002; Drake et al. 2002). 

Kaplan, van Kerkwijk, & Anderson (2002) reexamined the WFPC2 images, using a 
more sophisticated point spread function-fitting technique to measure source positions and 
improved (and unpublished) geometric distortion corrections. They concluded that the par- 
allax is 7±2 mas, resulting in a distance about twice what Walter (2001) measured. If true, 
the observational support for an extremely soft equation of state is removed. 

Three observations are the minimum required to measure a parallax. Since the expected 
parallax corresponds to a sub-pixel shift in the PC camera (45.5 mas pixel"^), we sought and 
were awarded a fourth observation, which occurred on 2001 March 24, in order to confirm 
the parallax. 

2. The Data 

The fourth WFPC2 image was scheduled near the time of maximum parallactic displace- 
ment on March 30. The first WFPC2 observation was described by Walter & Matthews 
(1997), and the next two were discussed by Walter (2001). The fourth observation is a 
7400 second observation consisting of 6 exposures at two dither positions with nominal off- 
sets of 5.5 pixels along each axis. The nominal pointing position was the same as in all of 
the other observations. The F606W filter was used to maximize the number of astromet- 
ric comparison stars. The choice of the same filter and pointing position was designed to 
minimize differential instrumental distortions. 

The particulars of all four observations are summarized in Table 1. The spectropho- 
tometry has been discussed by Pons et al. (2002). 



-3- 



3. Data Analysis 

We reanalyzed the images obtained at all four epochs, taking into account corrections 
discussed by Kaplan, van Kerkwijk, & Anderson (2002). We made no attempt to perform 
absolute astrometry. All positions are measured with respect to the first epoch (1996.7). 
The images were re-downloaded from the Multimission Archive at STScI (MAST) prior to 
analysis to ensure that the best instrumental cahbrations were applied using the MAST 
on-the-fly-calibration facility. 

We employed two independent measurement techniques and three independent analysis 
techniques, as described below. 

3.1. Measurement of Source Positions 

The first measurement of the source positions was performed by coadding and median- 
filtering the data obtained at each dither point, resulting in two images per epoch. We did 
not analyze the individual images, primarily because of S/N considerations. We fitted the 
positions of the targets in each pair of images, with a 2-dimensional Lorentzian function 
as template (using the IDL MPFIT2DPEAK^ function). We corrected the raw Y positions 
for the 34*'^ row error, using the prescription in Anderson & King (1999), and then apphed 
the geometric distortion correction from Holtzmann et al. (1995). Uncertainties in the po- 
sitions are the formal la uncertainties of the fit parameters, based on counting statistics 
in the images. To these uncertainties we added in quadrature an uncertainty of 0.03 pixels 
to account for systematic effects of non-uniformities in the intra-pixel response. For the 
subsequent analysis, we used the same reference objects (mostly field stars) as employed in 
Walter (2001), except that stars 109, 121, 125 and 126 were removed because of difficulties 
in fitting their positions. 

In the second measurement method we used the HSTphot software (Dolphin 2002). This 
code fits the point spread functions generated with the TinyTim (Krist 1993) software. We 
fitted all the images from each observation, at both dither positions, simultaneously. Dolphin 
(2002) claims an astrometric accuracy of 0.03 pixels with this software. The HSTphot 
astrometry corrects for the 34*'^ row error, but does not account for the geometric distortions, 
for which we applied the Holtzmann et al. (1995) correction. We followed all the processing 
steps described in the HSTphot manual prior to running HSTphot. HSTphot confirms that 
the objects identified as extended in Table 2 of Walter (2001) are indeed extended; all other 
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objects are consistent with point sources. We did not include the extended objects 104, 108, 
or 109 in this analysis. 

3.2. Analysis of the source positions 

We analyzed the measured positions by three independent methods. The derived proper 
motions and parallaxes are presented in Table 2. 

3.2.1. Full Astrometric Solution 

First, we performed full minimizations for the proper motions and parallaxes of the 
objects in the field, including image off^sets, residual rotations from the nominal roll angles, 
and scale factor changes from the nominal plate scale (45.5 mas pixel"^). This procedure 
was performed twice, both excluding and including the neutron star in the optimization. In 
the former case, the proper motion and parallax of the neutron star were obtained using 
the image offsets, residual rotations and scale factor changes from the other objects. This 
distinction was made since the expected proper motion and parallax of the neutron star are 
much greater than for those expected from the field objects. As anticipated, we found the 
parallax and proper motion of the neutron star are slightly smaller in the second analysis. 
The results quoted in Table 2 are from the first analysis. 

3.2.2. Independent X, Y Regressions 

Secondly, we determined the proper motion and parallax of the neutron star indepen- 
dently in the N-S and E-W directions. We registered the images with the assumption that 
the mean proper motions and parallaxes of the field stars are negligible. We rotated the 
measured positions to an equatorial coordinate frame using the nominal roll angles. We 
registered the images by shifting by the median offset in each coordinate. We iterated the 
registration, excluding stars whose residual differences are significant at >3(7 significance. 
Registration using a weighted mean shift produced insignificant differences. We then deter- 
mined the deviations from the nominal roll angle and plate scale by minimizing the diffences 
between the positions at each epoch and those of the first epoch. After resetting the roll an- 
gles and the plate scales, we re-registered the images. Uncertainties in the image registrations 
are about 0.02 pixels (1 mas) in each coordinate. 

We then shifted the measured position of the neutron star by the plate offsets, residual 
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rotations and plate scale changes. The parallax vector was determined by independent linear 
regression in both right ascension and declination. The true parallax is the projection of this 
vector in the direction of the parallactic motion (position angle 83°). 

In both the independent X,Y regressions and the full astrometric solution we find that 
the differences from the nominal rotation (< 0.02 degrees) and plate scale (< 0.03%) are 
small but significant. These affect the determination of the parallax at about the 20% level. 

3.2.3. Proper Motions From Annual Pairs 

As a safety check, we also determined the proper motions of all objects using the pairs 
of observations separated by integral years. The residuals at the half-year intervals are the 
sum of the parallactic shift and the measurement errors. We found a residual shifts in right 
ascension consistent with the parallax of the neutron star determined in the other analyses, 
but no significant residual shift in declination. The proper motions and parallaxes determined 
in this measurement are fully consistent with those in Table 2, but with uncertainties about 
a factor of three larger. 

4. Discussion and Astrophysical Implications 

The various techniques yield the same results within the uncertainties, and are presented 
in Table 2. For the subsequent discussion, we use the 117 ± 12 pc distance obtained with 
full least squares minimization using the HSTphot positions, which have slightly smaller 
nominal errors than those of the 2D fits. This result is robust because it is based on 4 
observations, two at each apex of the parallactic ellipse. This parallax is significantly smaller 
than that published earlier (Walter 2001). We attribute the difference to the omission of the 
geometric distortion of the camera, which amount to nearly 4 pixels at the edges of the fields 
(Holtzmann et al. 1995). The revised parallax is shghtly larger than, but agrees to within 
errors, with that reported by Kaplan, van Kerkwijk, & Anderson (2002). 

We examined all the stars in the field for their proper motions and parallaxes. The mean 
proper motions in (X,Y) are (-l,2)xl0~^ mas yr~^, respectively, and the mean parallax is 
4xl0~^ mas. This justifies the assumptions used for the image registration in §3.2.2. 
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4.1. Origin and Age 

(Walter 2001) suggested that RX J185635-3754 had its origin in the Upper Scorpius 
OB association. This conclusion is unaltered by the present analysis. The proper motion is 
essentially unchanged from what was originally reported, but the tangential space velocity 
is revised to 185(D/117 pc) km s~^. For assumed radial velocities larger than -150 km s~^, 
the projected path of the neutron star traverses the projected position of the Upper Scorpius 
OB association within the last 2 million years. 

Assuming that the neutron star originated in this association, the distance of closest 
approach of the neutron star to the center of the association is a function of the assumed 
radial velocity. For a present distance of 117 pc, the smallest separation of the star from the 
center of the association is about 8.5 ± 2.0 pc, which occurs for an assumed radial velocity 
of -10 km s~^. The closest approach occurred in this case about 0.5 million years ago. This 
separation is appreciably smaller than the size of the association, extrapolating its present 
size backwards in time. 

If RX J185635-3754 was born in the Upper Sco association, and its age is 5x10^ years, 
it is no longer a viable candidate for the binary companion of the runaway O star ( Oph. 
Hoogerwoerf, de Bruijnc, & de Zeeuw (2001) argue that the pulsar PSR J1932+1059 is the 
likely companion. RX J185635-3754 is either the result of a more recent supernova, or is 
unrelated to the Upper Sco association. 

The conclusion of Pons et al. (2002) that RX J185635-3754 is, within uncertainties, on 
the standard cooling curve expected for neutron stars is unaltered by the smaller revised 
age: the luminosity of the neutron star is also increased. The estimates quoted in Pons et 
al. (2002) for the magnetic field strength and spin period of the neutron star, based upon 
the observation of a bow shock nebula by van Kerkwijk & Kulkarni (2000) and a discussion 
by Walter & Wijers (2002), are revised to S ~ 10^^ G and P ~ 0.15 s. 

4.2. The Radius 

The most important reason to measure an accurate distance to this neutron star is to 
constrain its mass and radius. Our measurements indicate that, to a first approximation, 
these quantities will be approximately twice those estimated by Pons et al. (2002). For the 
heavy-element model atmosphere fits (Pons et al. 2002) of the spectral energy distribution 
of the neutron star, from optical to X-ray wavelengths, the revised mass and radius of 
RX J185635-3754 are shown in Figure 1. The model atmosphere fits, coupled with the 
revised distance of 117 pc, constrain both ~ 15 ± 3 km and the redshift z ~ 0.35 ± 0.15. 
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In turn, the quantities R and M are constrained to i? ~ 11.4 ± 2.0 km and M ~ 1.7 ± 0.4 

Mq. These values are permitted by a large number of current equations of state (Lattimer 
& Prakash 2001), including those containing exotic matter such as quarks in their cores. 
However, these constraints appear to be inconsistent with extremely soft equations of state. 

Since the analysis by Pons et al. (2002), new high resolution Chandra X-ray spectra have 
become available. Burwitz et al. (2001) find that the Chandra X-ray spectra is consistent with 
a blackbody temperature of 63 eV, and show that there is no evidence for the absorption 
lines and edges expected from a non-magnetized heavy-element atmosphere. The hotter 
temperature requires a smaller angular diameter R00/-D ~ 0.037 km pc^^ for the X-ray 
emitting blackbody (the statement by Drake et al. 2002 that the predicted optical flux 
increases by 10% is incorrect; the predicted optical flux actually decreases by about a factor 
of 3 relative to the extrapolation from the ROSAT blackbody because the inferred angular 
diameter decreases by nearly a factor of 2). This imphes it!oo — 4.2 km, even including a 
revised, increased distance (Drake et al. 2002). On this basis, together with the featureless 
spectrum and lack of significant pulsations, Drake et al. (2002) suggested that the object 
may be a self-bound quark star. 

However, as demonstrated by Pons et al. (2002), and earher by Pavlov et al. (1996) and 
by Rajagopal & Romani (1996), the inferred radius depends critically upon the details of 
the atmosphere and the spectral energy distribution. Comparing the ROSAT and EUVE 
X-ray observations with optical and ultraviolet data, Pons et al. (2002) showed that the 
optical and ultraviolet radiation cannot originate from the X-ray blackbody (see Figure 2). 
Either the X-ray emission arises in a hot polar cap, or there is significant modification of 
the spectral energy distribution from radiative transfer through a stellar atmosphere. The 
ROSAT/EUVE thermal component has an effective blackbody temperature of 55 eV, so 
the higher temperature inferred from Chandra observations only makes this argument more 
striking. 

If the surface has two different thermal components, we can follow the formalism of Pons 
et al. (2002) to estimate the true radius. Extrapolation of the 63 eV blackbody fit to the 
Chandra spectrum underpredicts the optical flux by a factor of 6 (Figure 2). Attribution of 
the optical flux to the cooler part of a two-temperature blackbody surface results in Rqo/D 
<0.18, or a maximum radiation radius Rqo ^ 21 km at a distance of 117 pc. In reality, the 
situation will be more complex than a simple one or two component blackbody. Nevertheless, 
the atmospheric models of Pons et al. (2002) have the net effect that the spectral emissions 
from a neutron star behave similarly to that of the two component blackbody model: the 
total radiating surface area must be substantially larger than a one component blackbody 
model to reconcile the optical data and X-ray data. On this basis, we suggest that the 
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Chandra data will not dramatically change the conclusions summarized in Figure 1. 

We are examining model atmosphere fits to the Chandra spectrum as well as to the full 
multiwavelength spectral energy distribution (work in progress). Pons et al. (2002) discussed 
only non-magnetic model atmospheres. New generations of magnetic model atmospheres 
are becoming available, and may be able to better reproduce the observed spectral energy 
distribution, as well as the lack of significant spectral features and pulsations. As we cannot 
yet select which atmospheric model best represents the full spectral energy distribution, the 
radius inferred from blackbody fits to the X-rays alone represent no more than lower limits 
to the true radiation radius. 

Wc emphasize in addition that the apparent disagreement between the Chandra and 
ROSAT spectral fits is troubling, as is a possible disagreement between the Chandra and 
EUVE fluxes where their responses overlap. Pons et al. (2002) found that the ROSAT and 
EUVE fluxes were mutually consistent, in comparison. There may be unresolved calibration 
issues which affect our ability to determine the radius at the level needed to distinguish 
various equations of state. Until we understand the full spectral energy distribution, it is 
premature to infer that the radius is substantially smaller than that expected from a normal 
neutron star. 
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081490197A from the Space Telescope Science Institute, by LTSA grant NAG57978, and by 
the USDOE grant AC02-ER-40317. We acknowledge stimulating converations with F. Ozel 
concerning magnetized neutron star atmospheres, and appreciate ongoing discussions with 
M. Prakash about the implications for nuclear theory. 
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Table 1. WFPC2 Observation Log 



Program 


Date 


Root 


Roll 


Duration 




ITT 




dog 


s 


6149 


1996 Oct 6 


U3IM01 


129.5 


4800 


7408 


1999 Mar 30 


U51G01 


-51.6 


7200 


7408 


1999 Sep 16 


U51G02 


124.2 


5191 


8567 


2001 Mar 25 


U62501 


-52.78 


7400 



Table 2. Astrometric Solutions 



Technique Proper Motion Position Angle Parallax Distance 

mas yr~^ deg mas pc 



HSTphot positions 

full matrix inversion 332.3 ± 0.4 100.45 ± 0.04 8.5 ± 0.9 117 ± 12 

independent X,Y regressions 332.7 ± 0.5 100.35 ± 0.05 8.8 ± 0.9 114 ± 12 

2D fit positions 



full matrix inversion 331.0 ± 0.7 100.47 ± 0.07 8.5 ± 1.5 117 ± 20 
independent X,Y regressions 331.0 ± 0.6 100.28 ± 0.08 9.1 ± 1.1 110 ± 13 
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Fig. 1. — Mass-radius diagrams for the uniform-temperature heavy element atmosphere 
models, revised from Figure 17 of Pons et al. (2002), for an assumed distance of 117 pc. 
Upper and lower panels are for Fe and Si-ash compositions, respectively. Solid and dashed 
curves are for equations of state labelled following Lattimer & Prakash (2001). The dashed 
line labelled "causality" is the compactness limit set by requiring equations of state to be 
causal. Dotted lines are contours of fixed i?oo- The crosses denote the masses and radii of 
models which best fit the optical and X-ray data at the indicated distance, and the hatched 
regions surrounding them include the nominal errors indicated in the constraint relations in 
equations (4)-(7) of Pons et al. (2002), as well as the nominal error in the distance. 

Fig. 2. — The spectral energy distribution of RX J 185635-3754. The Chandra LETG spec- 
trum is plotted along with the the EUVE and HST fiuxes from Pons et al. (2002). The dashed 
curve is the best blackbody fit to the Chandra LETG data from Burwitz et al. (2001); it 
underpredicts the optical fluxes. The dotted curve is the sum of the X-ray blackbody plus a 
15 eV blackbody with an angular diameter flve times larger than the X-ray blackbody. The 
two flts are indistinguishable at short wavelengths. 
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